Momentum-resolved tunneling between Luttinger liquids 
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We study tunneling between two nearby cleaved edge quantum wires in a perpendicular mag- 
netic field. Due to Coulomb forces between electrons, the wires form a strongly-interacting pair of 
Luttinger liquids. We calculate the low-temperature differential tunneling conductance, in which 
singular features map out the dispersion relations of the fractionalized quasiparticles of the system. 
The velocities of several such spin-charge separated excitations can be explicitly observed. Moreover, 
the proposed measurement directly demonstrates the splintering of the tunneling electrons into a 
multi-particle continuum of these quasiparticles, carrying separately charge from spin. A variety of 
corrections to the simple Luttinger model are also discussed. 

PACS; 71.10.Pm, 73.21.Hb, 73.40. Gk, 73.63.Nm 



The only universally accepted example of a non-Fermi 
liquid metallic state is the one dimensional (ID) Lut- 
tinger liquid (LL) Remarkably, the quasiparticle exci- 
tations of a LL are fractionalized, comprising a diverse set 
carrying spin separately from charge, and charge in frac- 
tions of the electron charge e. LL behavior has been ob- 
served experimentally in carbon nanotubes through 
strongly energy dependent local tunneling, and more re- 
cently in GaAs quantum wires through power-law res- 
onant tunneling lineshapes |^]. Evidence of charge frac- 
tionalization has also been seen in shot noise experiments 
using fractional quantum Hall edge states B, which are 
somewhat special chiral Luttinger liquids g . Despite 
these successes, no direct experimental evidence of frac- 
tionalization has ever been obtained in a non-chiral ID 
system. In this letter, we show that measurements of the 
non-linear tunneling conductance between parallel Lut- 
tinger liquids in a transverse magnetic field provide a 
direct spectroscopic probe of fractionalization. The re- 
sults, as described below, give very similar information 
to an ideal photoemission experiment. Indeed, some in- 
dications of spin-charge separation were seen in photoe- 
mission spectroscopy of the quasi-lD cuprate SrCu02 
Tunneling spectroscopy has, however, the advantages of 
being possible on a single, isolated ID system and with 
potentially much higher resolution than photoemission. 

Controlled tunneling experiments between two parallel 
wires have been recently conducted using cleavcd-edgc 
overgrowth by O.M. Auslacnder ct al j^. The experi- 
mental geometry we consider is indicated schematically 
in Fig. 1. The two "wires" are in fact confined surface 
states, and electrical contact is made only to the upper 
wire via a Two-Dimensional Electron Gas (2DEG). With 
L' ^ L, nearly the full electrochemical potential drop oc- 
curs between the shorter (left) segment of the upper wire 
and the lower wire. Because of the uniformity of the 
barrier, momentum along the wire is conserved during 
tunnehng B. 




FIG. 1. Schematic experimental geometry (from Ref. 7) 

We assume that the barrier is sufficiently high as to 
establish a quasi-equilibrium state on cither side of the 
barrier, treating tunneling across the barrier pcrturba- 
tively. A non-interacting model Hamiltonian for the sys- 
tem neglecting tunneling is then 
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where a = 1,2 labels the upper/lower wire, a ^tii l^-" 
bcls the electron spin, and Ua and fia are the electrostatic 
and chemical potential of the a*'' wire, respectively. We 
choose by convention to take /ia = in equilibrium, so 
that k-pa = y/^mUa is the Fermi momentum in wire a. 
Neglecting the energy dependence of the tunneling ampli- 
tude w, the Zeeman shift (see below), and a small energy 
shift due to orbital magnetic effects within each wire, the 
tunneling Hamiltonian in the presence of a magnetic field 
Bz (in the gauge Ay = Bx) is 
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Here the magnetic wavevector Q ~ l-nBdjcfiQ, d is the 
center-to-center distance of the wires, and 4>q = hc/e 



1 



is the flux quantum. The one-dimensional (ID) tunnel- 
ing current density, J = iew{c\^C2ae-^''^^ — ^2021^6"*'^'^) 
can be calculated directly from Fermi's golden rule. The 
result, whose gross features appear experimentally in 
Ref. 0, is shown in Fig. ||, taking ^1 = V, ^2 = 0. 
This diagram can be understood physically by consider- 
ing all processes by which an electron can be transferred 
between the two wires, moving from an occupied to an 
unoccupied state. Geometrically, Fig. 2 is obtained by 
drawing the locus of values of V, Q for which one of the 
four Fermi points lies on the other parabola when shifted 
vertically and horizontally by V and Q, respectively. 
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FIG. 2. Tunneling current in the V — Q plane for the 
non-interacting model. Regions of non-zero current are 
shaded, and dl /dV has delta-function singularities along the 
boundaries between shaded and unshaded regions. 

Zero bias features occur when (A) Q ~ ±(fcpi + k-p2), 
(B) Q = ±(A:fi — k^2)- It is in the low-energy region 
of the V — Q plane near these four features that univer- 
sal features arise in the interacting system, to which we 
now turn. Focusing on the the low-bias regime, we first 
decompose the electron fields into right and left movers, 
Caa{x) = Cflaa(a;)e''=''''^ + CLaa{x)e~'-^'^''^ , which are de- 
scribed by the Luttinger Hamiltonian (taking /^^ = for 
simplicity) 



Ha = 



' f) r 
Laa'-'x'^Laa 



(3) 



where Wpa = kj^a/fn- In general, these right- and left- 
moving Fermions undergo a diverse set of scattering pro- 
cesses mediated by the Coulomb interaction (screened 
by the 2DEG). A systematic study of these terms ||l^ re- 
veals an important simplification due to the experimental 
fact that the width W oi & typical cleaved edge quan- 
tum wire is large compared to the Fermi wavelength |^] 
[kj^aW » 1). We also expect the 2D screening length 
As ^ W . These properties imply a strong suppres- 
sion of two-electron backscattering processes. We there- 
fore adopt a forward scattering model retaining only the 
strongest (unsuppressed) interactions. 
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c[^aa'^Laa- H ncglccts momentum dependence of the 
forward-scattering interactions, and is hence valid for 
eV < hv-pa/^s- The interactions Vab can be roughly 
estimated as Vn ~ V22 ^ (2eVe) ln(2As/Ty), V12 ~ 
(2e^/e) ln(As/d) where e is the electron charge, e the wires 
dielectric constant. Experimentally, d — ^ |@], so 
the inter-wire interaction V12 is not negligible. Thus the 
tunneling conductance does not in fact probe the spec- 
tral properties of two decoupled IDEGs. Nevertheless, 
all properties of the interacting Hamiltonian Hq -\- iJjnt 
can be calculated exactly by bosonization (LL theory). 

From now on we will focus on the current at zero 
temperature around the point A of Fig. |^. An anal- 
ogous description can be obtained around each of the 
low bias points of fig. || [|lO|. To lowest order in per- 
turbation in w, we can write the current density J as 
the sum J = 2e\w\'^{J+ — J_), where J+^J- are pos- 
itive functions which satisfy (time-reversal symmetry) 
J+{Q,V) = J_(— Q,— F), and J+ (J_) is nonzero only 
for 1/ > (F < 0). Around point A, q ^ Q - (kp^ + kp^) 
is small, and 
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where C^^^ is the Euclidean correlation function 
C^^^{x,t) = (|c^iC^2(2;,T)cLcm(0,0)y LL theory 
then gives 
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where Ca = (^1)°^^- and ao (of order W) is a small dis- 
tance cut-off for the LL description. Fractionalization is 
evident formally in Eq. ^ through the multiple branch 
points characterized by distinct charge and spin veloci- 
ties, Wci/2 and Wsi/2- For two independent wires (V12 = 0), 
it is well-known that 77 = 1/4 and 6*1, 02 > 1/4. In the 
present case, the strong interactions between the wires 
make 77 interaction-dependent, and Oi, 02 can take values 
smaller than 1/4. Indeed in the case of two coupled iden- 
tical wires {v-pi = Wf2 = Vp^Vu = V22), these exponents 
take the values 9i 

Vl + 2Ct^i] 



1/(4^1 -f 2(1/11 +yi2)/(7n;F)), 
0. 



72 ~ Y i -I- zi^i/11 — Vi2)/(7rwF)/4, ry = 0. General but 
non-illuminating formula are postponed to . 

Letting x tu in Eq. ^ and integrating over t gives 
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where the complex function h{u) is defined by 

h{u) = [5- i{V + gu)]2«i+2»="iC7|-*^(w, i + e). 
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FIG. 3. Integration contour and cuts for V,q > 0. The 
original integral of h{u) defined on — oo, +00 (contour 1) is 
modified into a integral on contour 2+3. 

The contour of integration in (^) can then be deformed 
according to Fig. |[ As h{u) vanishes faster than at 
infinity, we are reduced to the integral of h{u) on the con- 
tour 3 (see Fig. ||) . More explicitly in the situation where 
I'si < Vci < i's2 < i'c2, and for e.g q > 0, the remaining 
integral can be written as the sum of real integrals: 



behavior is expected to be manifested for eV < ep 
10— 20ml^ in cleaved edge samples. Thus for experiments 
with T in the Kelvin range, the thermal rounding is min- 
imal except at very low bias. Second, for k^aW ^ 1, the 
dominant impurity process is elastic forward scattering. 
This is easily included by convolving G{V, q) (in q) calcu- 
lated above with a Lorentzian of half- width Aq = 1 / Hpi , 
and has several consequences. The singularities at > 
are rounded over an energy width I/tci ~ Vp^i/tei- Non- 
vanishing weight also appears outside the kinematically 
allowed region. At low bias and temperature, this weight 
is itself singular and displays a pronounced low bias con- 
ductance dip. In particular, for eV,kBT,v*q <^ ^/tci, 
G{V,q,£ei) - \w\^[iei/{l + qHli)]TPFp{eV/kBT)., with 
< /3 = 26*1 + 202 - K 1, where Fi3{X) is the well-known 
scaling function for point-contact tunneling, satisfying 
F^(0) = 1 and Fp{X) - Xf^ for X > 1. 



J+(Q, V) oc 2 sin(2^(0i + 02))r(l - 29^ - 262) 
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x^e{V ~ qva) sin <j>a \h(u)\du, (7) 

with vi = Vci,V2 = Vs2,vt, = Vc2,V4 = +00 and 0i = 
7r(6'i-?7),02 = 7r(6ii-?/+i),03 ^ 71(61+62 + 1/2). Similar 
expressions for J+ and J_ can be derived for any other 
order of the velocities, and for the other F ~ points of 
Fig.^ ijl^. A typical result is presented in Fig. ^. 

The density plot of the differential conductance (per 
unit length) G = dJ/dV in Fig. ^ directly exhibits ev- 
idence of electron fractionalization. First, non-analytic 
features appear along rays (3 per quadrant) whose 
slope gives the two charge and spin velocities. Sec- 
ond, dJ/dV is non-zero for any bias above the thresh- 
old \V\ > v*\q\, with V* = Min{vs2,Vci,Vc2)&{qV) + 
Min{vsi,Vci,Vc2)'d{—qV). Both these properties can be 
understood from kinematics. Each tunneling event cor- 
responds to a transfer of charge ±e and spin ±1/2 from 
one wire to the other, accompanied by the addition of mo- 
mentum q and energy eV. This creates a combination of 
(6 total) fractional "chargons" and "spinous" . The final 
state with the least energy for a given momentum q gives 
all the momentum to the slowest appropriate "particle" , 
which determines the threshold. Moreover, since the fi- 
nal states are six-particle excitations, kinematics allows 
any energy above v*]'?!, explaining the non-zero weight in 
dJ/dV as a multi-particle continuum. Detailed expres- 
sions for the singularities along the various rays in Fig. ^ 
can be extracted from Eq. |7| . 

We now turn to a discussion of the numerous effects 
left out of the above treatment. The two most signifi- 
cant corrections can be treated exactly. First, the spec- 
tral density is rounded on a scale set by temperature (the 
rounding can be calculated exactly to describe detailed 
lineshapes pG|). Experimentally, for e^/vp of 0(1), LL 
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FIG. 4. Top: Density plot of the ideal differential conduc- 
tance at point A of Fig. ^. Bottom: Differential conductance 
G{V, q — 1) versus potential difference V. Both plots are 
shown for values approximately appropriate to cleaved edge 
overgrowth: 61 = 0.2,02 = 0.32, r; — 0.18, with velocities 
Vsi = 1, Vci = 1.36, Vs2 = 2, Vc2 = 2.56 in units of dfi. 

Even at r = 0, effects not included in the LL model 
can broaden the spectral function. First, consider the ef- 
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fects of non-forward-scattering interactions between elec- 
trons at the edge. In the generic situation, fcpi 7^ kp2, 
and charge/spin-density- wave couphng (e.g. Hcdw = 

VcDW'0lji'0Li'0L2''^_R2) between the two wires is forbid- 
den. The dominant residual term is then the "ex- 
change" interaction (we reserve the usual term "backscat- 
tering" for impurity effects discussed below) within a 

single wire, e.g. Hox,i = -A^WpiV'L'^V'ffi ' fph^i'u^ 
where Ai characterizes the dimensionless backscattering 
strength in wire i. In the experimentally relevant situ- 
ation, Xi ~ (kpiW)'-^ ^ 1, so this is a weak interac- 
tion. Formally, such exchange interactions are known 
to be marginally irrelevant in the renormalization group 
sense, and a simple perturbative (bosonized "phonon" 
self-energy) estimate suggests a resulting "lifetime" scal- 
ing approximately linearly with energy u, (1/tox,i;)^ ~ 
Afw^/ (ci-|-C2Af ln^(eF,i/w)), where Ci/2 are order one con- 
stants. The magnitude of 1/Tex,i, however, is small, due 
both to the smallness of Xi and the additional logarith- 
mic suppression at low energies. Other "internal" cor- 
rections to the forward-scattering model, such as band 
curvature, are strongly irrelevant, and a similar self- 
energy estimate shows that they contribute only negli- 
gible 1/Tcurve ~ 0{uj'^ /ep) tcrms to the lifetime. 

Electrons in the wires also interact with those in the 
bulk 2DEG. Expressing the Coulomb interaction in the 
basis of wire and bulk states, one finds several distinct 
scattering channels. Most interesting are charge ±e de- 
cay processes, in which an electron/hole in one of the 
wires decays into an electron/hole in the 2DEG and an 
electron/hole pair. Because the 2DEG is a Fermi-liquid, 
however, strong phase-space restrictions reduce the rate 
of such decays at low energies: crudely, I/t^ ^ CeW^/ep, 
where the 0{1) constant Ce depends upon the density 



of the 2DEG, etc. Interestingly, if /cp 



< 3A;f i, such 



processes are kinematically forbidden (ce = 0). Simi- 
lar, albeit slightly less restrictive kinematic conditions 
reduce the phase space for charge ±2e decay, in which 
Cooper pairs (hole pairs) are scattered between the wires 
and 2DEG. All such processes in which charge is trans- 
ferred are further reduced by LL orthogonality catastro- 
phe effects. Finally, there are forward-scattering pro- 
cesses in which no charge but only energy and small mo- 
menta are transferred between the wires and 2DEG. Be- 
cause the Coulomb potential suppresses long-wavelength 
charge fluctuations in the 2DEG, we expect this also to 
be a weak effect, and will explore it in greater detail in a 
future publication |10|. 

In addition to the elastic forward scattering discussed 
above, disorder also mixes surface and extended states, 
allowing conduction at F = T = 0. A lower bound 
on the broadening of spectral features in the tunneling 
conductance is provided by the experimental determina- 
tion of the 1D-2D scattering length , iiD~2D ~ 6/i in 
Rcf. H; hence 1/tid-2D > t^F,i/^i_D-2D- The latter is a 



lower-bound to the elastic contribution, since there is also 
a small-momentum elastic component which contributes 
to the broadening but not Iid-2D (or the longer elastic 
"backscattering" length £b defined in Ref. ^). 

Two additional physical effects have not been taken 
into account above. First, Zeeman coupling leads exactly 
lllOi to two superimposed copies of the LL results with Q 

^±^iv-^+v-^) for 



shifted left and right by AQ^^^ 



point A and AQy^ 



= ±i 



2 -(Wpi - Wp/) for point B 
Second, direct 2DEG-wire tunneling provides a second, 
albeit weaker, conduction channel but can be dis- 
tinguished experimentally since it contributes at low-bias 
for the full range /i;F2 — fcp.buik < |Q| < fcF2 + fcF,buik instead 
of at the singular points A, B. 

Finally, we comment on the validity of the perturba- 
tive treatment of tunneling. For q 7^ 0, the tunneling 
perturbation is formally an (infinitely!) irrelevant opera- 
tor. Thus for sufficiently weak tunneling at non-zero q the 
perturbative treatment is accurate. At q = 0, however, 
tunneling is relevant: at an energy e, the effective tunnel- 
ing amplitude is w{t) = w(eo/e)^/^~^^~^^, where the cut- 
off eo ^ J. Equating this to cqj we estimate the cut-off 
energy eV* - eo(Weo)^/^^"^^'"^^''- For el/, Up,.g < eV\ 
the perturbative treatment breaks down. Physically, in 
this regime we expect coherent motion between the two 
wires. Due to the necessary complication of including dis- 
sipation in the leads, we therefore postpone the analysis 
of this limit to a future publication . 
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